Multimode quantum limits to the linewidth of an atom laser 
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The linewidth of an atom laser can be limited by excitation of higher energy modes in the source 
Bose- Einstein condensate, energy shifts in that condensate due to the atomic interactions, or phase 
diffusion of the lasing mode due to those interactions. The first two are effects that can be described 
with a semiclassical model, and have been studied in detail for both pumped and unpumped atom 
lasers. The third is a purely quantum statistical effect, and has been studied only in zero dimensional 
models. We examine an unpumped atom laser in one dimension using a quantum field theory using 
stochastic methods based on the truncated Wigner approach. This allows spatial and statistical 
effects to be examined simultaneously, and the linewidth limit for unpumped atom lasers is quantified 
in various limits. 

PACS numbers: 03.75.Pp, 05.10.Gg, 03.70.+k 



I. INTRODUCTION 

The experimental realization of Bose-Einstein conden- 
sates (BECs) provided a testbed for many fundamental 
questions in interacting quantum systems, as well as pro- 
viding a general tool for investigating aspects of atomic 
physics 0, 0, H, H| . One major application that BECs 
offer is the possibility of creating an atom laser. Atom 
lasers are the matter wave analogy of the optical laser, 
and are created by coupling atoms out of a BEC by us- 
ing some external means to change the state of a subset 
of the atoms in the condensate from a trapped to an 
untrapped state H, H, 0] . This can produce a beam of 
atoms that exhibits both spatial and temporal coherence 
B H EE HO, El EH- While the optical laser is well- 
studied and well-understood, the atom laser is still in its 
infancy, both in terms of a full theoretical description and 
experimental realization. 

An important property of a laser is its linewidth, which 
is of crucial importance for many measurements, includ- 
ing spectroscopy and interferometry [l4l |. For an "un- 
pumped" optical laser, which is simply a cavity contain- 
ing light, the linewidth is simply the inverse damping 
rate. For a pumped optical laser, the linewidth is often 
limited by technical effects of the pumping, but the fun- 
damental quantum limit is the Schawlow-Townes limit, 
which can be seen as a phase diffusion process arising 
from the addition of spontaneously emitted photons with 
random phases into the radiation field [151 ]. The deriva- 
tion of this limit in the optical case relies on the fact 
that photons do not interact with each other, and conse- 
quently the damping rate of the cavity is independent of 
the photon population. 

For an atom laser, things are more complicated. The 
different dispersion relations of free atoms compared to 
free photons means that the relationship between outcou- 
pling and linewidth of a cavity is non-trivial, except in 
the weak outcoupling limit • When the Bose-Einstein 
condensate is stable in a single lasing mode with no inter- 
actions, the equivalent of the Schawlow-Townes limit also 



exists for weak outcoupling 16]. Unlike photons, how- 
ever, atoms interact with each other strongly, resulting 
in a nonlinear interaction term arising from atom-atom 
collisions. This has several important affects on the dy- 
namics of an atom laser. Semiclassical calculations show 
that the single-mode operation of an atom laser tend to 
occur in parameter regimes where the interactions are 
dominant , EH . Interactions also cause diffusion of the 
phase of the condensate, and this quantum noise is the 
dominant contribution to the linewidth when the single 
mode approximation is valid [19]. In practice, this phase 
diffusion will also be spatially dependent, and the spatial 
effects may interact with the effects of the quantum noise 
on the output. 

This paper examines the linewidth of an unpumped 
atom laser with strong interactions, without making sin- 
gle mode or semiclassical approximations. The results 
are compared to simulations based on the semiclassical 
approximation and theoretical results based on a single- 
mode, zero-dimensional model. We find that our results 
scale similarly to the semiclassical and zero-dimensional 
models, but the quantitative details depend on the sys- 
tem parameters. 

General quantum field theoretic problems even in a sin- 
gle dimension are intractable by brute force methods, and 
the dominant interactions in atom laser systems means 
that linearized analytic calculations are not valid. An al- 
ternative to brute force calculation is offered by stochas- 
tic phase space methods. In this paper, we apply a trun- 
cated Wigner approach to the problem of determining 
the linewidth of an experimentally realistic, multimode 
atom laser. 



II. THEORETICAL METHODS 

Our model for the atom laser consists of a population 
of three level atoms in a trap, as shown in Figure [1] 
Bose condensed atoms in state |1) are confined by a har- 
monic trapping potential. These trapped atoms are then 
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coupled to an untrapped state |2) in which they do not 
see the confining potential and consequently leave the 
trap forming the atom laser beam. The coupling from 
the trapped to the untrapped state is done via a Raman 
transition using the intermediate level |3). Raman out- 
coupling has been used by two atom laser experiments 
@; H3| ! and leads to an atom laser with superior proper- 
ties such as a higher flux and higher brightness [2CJ], the 
ability to give the beam directionality [3, |2(|, and the 
possibility of creating non-classical states of the beam 
[2l| . The Raman outcoupling scheme reduces to the rf 
outcoupling scheme in the limit of a zero momentum kick. 
The three level system shown in Figure [1] can be reduced 
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FIG. 1: Atomic level scheme. Using a Raman process, 
trapped atoms in the condensate (|1)) are transferred to the 
untrapped state (|2)) via an intermediate state (|3)). 



to a two level system by adiabatically eliminating the 
upper level. The effective second-quantized Hamiltonian 
describing this reduced two-level system is given by 
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where *i(r) and 1 & 2 (r) describe the trapped and un- 
trapped matter fields respectively, S = A23 — A13 is the 
two photon detuning, Q = f^ 3 i7i2/Ai3 is the two-photon 
Raman Rabi frequency, and ko = k2 — ki is the momen- 
tum kick imparted to the untrapped atoms from the two 
photon transition. The trap has been assumed to be 
isotropic and the trapping frequency is u>. As we have 
included position dependence in the matter fields, the ef- 
fective Hamiltonian describes the full multimode nature 
of the problem, and also includes non-Markovian effects. 
We assume that the atomic gas is sufficiently cold and 



dilute, so only binary collisions are relevant and the non- 
linear potentials are defined by 



Uij = AnTi 2 aij /m, 



(2) 



where ay is the s-wave scattering length between atoms 
in state |i) and state 

When the quantum statistics have no effect on the dy- 
namics of the mean field, these equations can be solved 
semiclassically, using the Gross-Pitaevskii (GP) equation 
to describe the atoms. The coupled GP equations arising 
from Eq. JT]) are given by 
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where we have moved to a rotating frame. These equa- 
tions can be solved numerically in one, two or even three 
dimensions, depending on the spatial resolution required 
and computational resources available. 

In the absence of interactions, the field theory solution 
and the semiclassical solution are identical, and demon- 
strate that in the limit of weak interactions the output 
approaches a Lorentzian with a linewidth that depends 
on the dispersion relations of the output field [l4| . In this 
Markovian limit the linewidth is the inverse of the time 
taken to empty the condensate. While the proportional- 
ity constant changes depending on the potential seen by 
the output field, in the limit the linewidth is always pro- 
portional to the rate of the state-changing mechanism, 
and therefore has no lower bound. At higher coupling 
rates, the timescale of the transport of the atoms from 
the outcoupling regime becomes comparable to the back- 
coupling rate, and the beam demonstrates increasingly 
complicated spatial and spectral behavior, up to and in- 
cluding the shut down of the outcoupling process alto- 
gether [22I ]. 

In the presence of interactions, the semiclassical equa- 
tions given by Eq. (J3]) can be solved numerically, and 
demonstrate important issues even in the limit of weak 
outcoupling, such as a 'chirp' in the frequency of the out- 
put beam due to the decay of the energy of the trapped 
state PI . The 'chirp' can be removed by adjusting the 
two-photon detuning during the outcoupling, but the 
presence of interactions also cause the quantum field to 
exhibit effects not seen by the semiclassical model. The 
phase diffusion of the atomic field will cause a broaden- 
ing of the linewidth of the atom laser, and in the limit of 
weak outcoupling, this will be the dominant lower bound 
to the linewidth. It is also possible for these effects to 
have spatial effects in the high outcoupling limit, as the 
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three timescales in the problem become comparable to 
each other. 

We will first estimate the size of the effects of the quan- 
tum statistical effects before calculating them in a mul- 
timode model. 



A. Single mode limit 

Weakly outcoupled atom lasers use condensates that 
should be well described by the Thomas- Fermi limit. In 
this limit, the chemical potential depends on the dimen- 
sion of the condensate, and is given by 
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where the nonlinear potential U is given by Eq. ([2]), u> is 
the geometric mean of the trapping frequencies, and A 
and I are dimensional reduction factors corresponding to 
an area and a length respectively. 

These dimension reduction factors arise when one has a 
condensate that is tightly confined in one or two dimen- 
sions, corresponding to a highly anisotropic trap. The 
tight confinement along one or two axes ensures the non- 
linearity has negligible effect in these directions, meaning 
the matter field operator will factorize with the trans- 
verse dependence completely described by a coherent 
state occupation of the lowest radial trap mode. This as- 
sumption leads to a second quantized Hamiltonian iden- 
tical to Eq. (p}, except that the matter fields are one- or 
two-dimensional, the integral is over a length or an area 
rather than a volume, and the nonlinear potentials are 
scaled by a transverse area or transverse length. This 
area or length corresponds to the cross sectional area or 
width of the condensate along the tight trap directions 
[23l ] . Dimensional reduction is a useful technique for nu- 
merical simulation, as reducing a three-dimensional prob- 
lem to a one-dimensional problem greatly reduces the 
computation required. 

If one assumes that the condensate is in a coherent 
state, it has one standard deviation number uncertainty 
given by ^/N. Consequently the energy uncertainty in 
the condensate due to number fluctuations will be given 
by 
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which means the energy uncertainty of the condensate is 

(8) 
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Equations (|5])-(fT0]) give a rough measure of the funda- 
mental linewidth of an atom laser, provided that the 
linewidth limit is due to number fluctuations in the con- 
densate being transformed into phase fluctuations due to 
the nonlinear atom-atom interactions. In order to test 
the accuracy of this estimate, we require a fully multi- 
mode model of the system, for which we will use stochas- 
tic methods. 



B. Stochastic methods 

The standard approach to using stochastic methods for 
a specific problem is to express the density matrix of the 
system in a specific basis as a quasiprobability distribu- 
tion such as the P, positive P or Wigner distribution. 
The master equation of the system is then converted to 
an equivalent partial differential equation in terms of the 
distribution chosen, which can be cast in form of the 
Fokker-Planck equation 
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where a is a vector of complex fields and D = BB T is 
the diffusion matrix. Given a Fokker-Planck equation of 
the form (jlip . a fully equivalent formulation is given by 
the system of Ito stochastic equations 



-^- = A i + B ji r ]i (t) 



(12) 



where the T)i(t) are set of noise sources with zero mean 
that are delta-correlated in time. It is therefore possible 
to describe the evolution of a quantum system in terms of 
a set of Langevin or stochastic partial differential equa- 
tions. These stochastic equations contain all the informa- 
tion of the original quantum problem in that various en- 
semble averages of the stochastic variables correspond to 
various expectation values of the original quantum fields. 
In the case of the Wigner distribution, the stochastic av- 
erage of an expression consisting of complex- valued fields 
corresponds to the expectation value of the symmetrized 
version of the quantum field operator version of that ex- 
pression. 
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This means that if we wish to find the expectation 
value of some normally ordered operator O given by 



6 (a, a*) = c nm a) r 



(13) 



we need only calculate the stochastic average of some c- 
number function O s (a* , a) . For the Wigner distribution, 
O s is given by [l5[ 
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where the overbar indicates a stochastic average. Thus, 
for example, the stochastic average equivalent to the 
number operator is given by 



(fit a) = a*a . 



(16) 



The master equation corresponding to the effective 
Hamiltonian |T]) is given by 
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We determine this master equation in terms of the field 
operators Wi(r) and ^ 2 (r), and make the following re- 
placements to transform it into an equation of motion for 
the functional Wigner distribution: 
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ables tpi and i/> 2 . This gives 
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where K. = — h 2 V 2 /2m, V = mu> 2 r 2 /2 and where the 
third order functional derivatives have been dropped. 
These third order terms do not have a simple mapping 
to stochastic partial differential equations, and can be 
assumed to be negligible when the field has a high oc- 
cupation number. This "truncated" Wigner approach 
has been used successfully in a range of calculations 

[H [M m, m, H3, M M, m ■ 

A Fokker-Plank equation can be sampled with equiv- 
alent stochastic equations where the noise terms depend 
only on the second order derivatives. As Eq. (|23p has 
no second order terms, the equivalent equations are com- 
pletely deterministic once a particular initial state is cho- 
sen. Although the evolution is completely deterministic, 
we have not removed all effects of the quantum noise, as 
the noise will still enter in the choice of initial states. 

Comparing Eqs. (fTT|). ([12"]) and we see that Eq. 
([25)1 is equivalent to the following pair of partial differ- 
ential equations 
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We then rearrange the result in the form of a Fokker- 
Planck equation in terms of the c-number stochastic vari- 



where AV is the volume element of the discretization 
of the problem. The terms inversely proportional to 
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these volume elements compensate for the mean field of 
the vacuum, which is non-zero in the Wigner approach. 
Apart from these terms, these equations are identical 
to the semiclassical equations for the system given in 
Eq. ©. 

To determine the noise distribution we need to apply 
to our initial states, we make the assumption that the 
condensate is in a multimode coherent state, which is 
equivalent to assuming that each of the single mode fields 
at each of the grid points Xi is in a single mode coher- 
ent state \ipj(xi)). We can find an appropriate ground 
coherent state for the trapped atoms by integrating the 
semiclassical equations in imaginary time, and the out- 
put is assumed to be in the vacuum state initially. 

The Wigner distribution for a single mode coherent 
state | a) is given by 

W{yi,y 2 ) = ^ exp [-2 ((a* - a r f + (x 2 - a,) 2 )] (26) 

where a r and a.i are the real and imaginary parts of the 
coherent amplitude a, and y\ and y 2 correspond to the 
real and imaginary parts of the atomic field. This Wigner 
function represents a Gaussian uncertainty with a stan- 
dard deviation of one half along both the y\ and y 2 axes. 
Thus, to give our initial fields the correct statistics, at 
each grid point we need to apply noise of the following 
form: 



(27) 
(28) 



where the rjj(xi) are complex Gaussian noise functions 
with a standard deviation in the real and imaginary com- 
ponents of one half. 

To investigate the linewidth behavior, we solved the 
stochastic equations (|24l) and (f25l) numerically, using the 
mathematical package XMDS [3l|]. This open source 
package allows the solution of systems of both stochastic 
and deterministic partial differential equations on grids 
of arbitrary dimension. It also allows the use of MPI 
methods to distribute the computation over many pro- 
cessors, an approach that is crucial when carrying out 
multi-trajectory simulations such as those required by 
stochastic PDEs. These stochastic solutions were then 
compared with the semiclassical solutions obtained by 
solving Eq. 

The linewidth of the the atom laser beam is the width 
of the momentum distribution of the beam as a function 
of time. Using Eq. (JTHJ) along with the linearity of the 
Fourier transform, one can show that the expression for 
the beam momentum in terms of the stochastic variables 
is given by 
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The system was mostly solved in one dimension for 
reasons of computational efficiency, although a limited 
number of two-dimensional simulations were carried out 
in order to ensure that the scaling laws carried over from 
one dimension to two as expected. The dimensional re- 
duction of the problem from three dimensions to one was 
performed as described in Section Ul Al We chose a trans- 
verse area of A = 1.2 x lCU n m 2 for reduction to a ID 
problem and a transverse length of I = 3.46 x 10 -6 m for 
the reduction to a 2D problem. 

We focus on the weak outcoupling regime in these sim- 
ulations so that we can examine the effects of interaction- 
induced phase diffusion in isolation from other limits to 
the linewidth. In the strong outcoupling regime a num- 
ber of effects can cause density fluctuations in the con- 
densate which are reproduced in the output beam [22] |. 
and this classical noise dominates the resulting linewidth. 
Also, the chemical potential decreases with the number 
of atoms in the condensate, resulting in a smaller mean 
field kick to the atoms as they leave the trap. This causes 
a "chirp" in the momentum distribution of the beam, ar- 
tificially broadening the linewidth [14j . In the weak out- 
coupling limit, the interaction-induced phase diffusion is 
the dominant limit to the linewidth. 



III. RESULTS 

The quantum noise has no effect on the short term 
behavior of the atom laser. The main short-term fea- 
tures are shown in Figure [H which shows the outcoupled 
field in momentum space. In the first few microseconds, 
the output momentum spectrum mimics the trapped mo- 
mentum spectrum shifted to be centered around the total 
momentum kick from the Raman transition hko- As the 
atoms are accelerated out of the condensate due to the 
mean field repulsion, this peak moves, and over time the 
energy width of the output field reduces and this peak 
narrows. Thus, just after outcoupling begins, there is a 
spread of momenta in the beam ranging from the Raman 
kick hko up to the Raman kick plus the momentum corre- 
sponding to the energy fj, gained from the mean field. As 
time progresses, many more atoms accumulate in the part 
of the beam that is outside the condensate rather than 
in the outcoupling region within the condensate. Conse- 
quently, in momentum space, most of the momentum is 

now concentrated in a peak at Tik -- 



where A 14 is the volume element of the discretized grid 
in momentum space. 



h /cq + 2m/z, and 
it is the width of this peak that we will consider to be 
the linewidth of the atom laser. Early snapshots of this 
behavior are shown in Figure [2l 

Previous semiclassical analyses of the width of this 
peak have shown that it is essentially Fourier limited: 
that is, its width is inversely proportional to the out- 
coupling time 14| . This narrowing continues indefinitely 
according to a semiclassical analysis, a result that must 
inevitably fail at some point due to quantum statistical 
fluctuations of the phase of the BEC providing the atom 
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FIG. 2: Momentum space power spectrum of the output beam 
after 5ms (left) and 10ms (right). Over time, the atoms 
build up outside the coupling region, and gain a small ex- 
tra momentum kick due to the BEC-beam repulsion. Pa- 
rameters: N = 5 x 10 6 , uj = 250rads _1 , a = 1 x 10~ 9 m, 
fco = 2 x 10 7 m _1 . 



laser linewidth using both a semiclassical and a stochas- 
tic approach, we notice a striking difference in long-term 
behavior. The linewidth of the beam over time is plot- 
ted for two different parameter regimes in Figure U and 
Figure [5] These Figures show that both the semiclassical 
and stochastic results agree initially on the rate at which 
the linewidth narrows, but at later times the stochastic 
simulation shows that the linewidth hits a limit, while in 
the semiclassical simulation it continues to narrow. 

In the long time limit the shape of the semiclassical 
(dashed) curve in Figures H] and [5] is linear with a slope 
of -1, indicating that linewidth is inversely proportional 
to the outcoupling time. This agrees with the Fourier 
arguments in The fundamental linewidth limit due 
to interaction-induced phase diffusion for an approximate 
single mode model, given by twice Eq. ((5|), is shown on 
the figures as a horizontal bar, and agrees closely with 
the results of the multimode stochastic simulations. 



source. The quantum statistics become important on a 
timescale depending on the strength of the interactions, 
which was typically on the order of tens of milliseconds 
in our simulations. 

The momentum space power spectrum of the out- 
put beam was calculated by using Eq. (j2"5)) and taking 
the stochastic average over 1024 paths unless otherwise 
stated. The resulting power spectrum of the beam was 
then fit to a Gaussian envelope function and the full 
width, 1/e height taken as the linewidth. This defini- 
tion of linewidth corresponds to twice the standard de- 
viation of the Gaussian, and as such is given by twice 
the energy uncertainty given in Eqs. (f5|)- (fTTj)) . Due to 
the finite number of paths and the stochastic nature of 
the problem, the Gaussian fit is not exact, but it still 
resembles a Gaussian very closely. An example is shown 
in Figure [3J When we examine the narrowing of the 
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FIG. 3: Power spectrum of the output beam at t = 0.35 s 
according to the stochastic simulation (solid line) as well as a 
Gaussian fitting profile (dashed line). Parameters: N = 10 7 , 
uj = 250rads"\ a = 3 x 10" 9 m, k = 10 7 m- 1 . 
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FIG. 4: Linewidth narrowing as a function of outcoupling 
time. Dashed line shows the semiclassical result; solid line 
shows the result of the stochastic simulation. Horizontal line 
indicates the fundamental linewidth limit according to twice 



Eq. J8j). Parameters: N 
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To check the possibility that the good agreement be- 
tween the linewidth limit given by ([H]) and the results of 
the stochastic simulation is related to the restricted di- 
mensionality of the one dimensional simulations, we re- 
peated the simulation in two dimensions using 256 tra- 
jectories. In two dimensions the momentum distribution 
is now no longer a simple peak, but rather an arc in mo- 
mentum space satisfying the energy conservation relation 
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(30) 



where x is the transverse degree of freedom that we have 
added. An example of this two-dimensional momentum 
space density is shown in Figure [6l The peaks in the mo- 
mentum distribution are due to transverse structure aris- 
ing from interference effects in the beam and the outcou- 
pling process. The parabolic arc has a specific thickness 
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FIG. 5: Linewidth narrowing as a function of outcoupling 
time. Dashed line shows the semiclassical result; solid line 
shows the result of the stochastic simulation. Horizontal line 
indicates the fundamental linewidth limit according to twice 
Eq. ©. Parameters: N = 10 7 , w = 250rads _1 , a = 3 x 
10~ 9 m, fc = 10 7 m _1 . 



FIG. 7: Results for a two dimensional simulation. Linewidth 
narrowing as a function of outcoupling time, with energy 
spread measured in the longitudinal direction. Dashed line 
shows the semiclassical result; solid line shows the result of 
the stochastic simulation. Horizontal line indicates the fun- 
damental linewidth limit according to (|9}. Parameters: N = 
5 x 10 6 , u = 2000 rads -1 , a = 1 x 10~ 8 m, k = 2 x 10 7 m _1 . 



that narrows, and corresponds to the energy linewidth 
of the output beam. The results of this narrowing as a 
function of time are shown in Figure [7J Again, the fun- 
damental lower linewidth limit closely tracks that pre- 
dicted by @. Consequently the argument that number 
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FIG. 6: Momentum space density of a two dimensional atom 
laser, with p x = hk x , p z = hk z . The parabolic shape is due 
to energy conservation as per Eq. (|30[) . with the thickness of 
the line tracing out the parabola giving the linewidth. Pa- 
rameters: N = 5 x 10 6 , uj = 1500 rads -1 , a = 1 x 10~ 8 m, 
fc = 4 x 10 6 m _1 . 



fluctuations coupling to energy fluctuations cause a fun- 
damental limit to the linewidth appears to scale correctly 
across dimensions, indicating that Eqs. ([8])- (fT0|) will cor- 
rectly estimate the fundamental linewidth of a real atom 
laser experiment in three dimensions. 



IV. NUMBER SQUEEZING 

As the linewidth of the atom laser scales with the size 
of the number uncertainty of the condensate, it is natural 
to consider to what extent the linewidth could be reduced 
by minimizing number uncertainty via number squeezing. 

The most straightforward approach is to consider 
quadrature squeezing. In the single-mode case one de- 
fines amplitude and phase quadrature operators by 



(31) 
(32) 



where 4> is the phase angle at which the measurement is 
carried out. The variances of are unity for a coherent 
state, and consequently a state is squeezed if the variance 
of one the quadrature operators is less than one. Squeez- 
ing effectively repartitions the unavoidable uncertainties 
associated with simultaneous measurement of a pair of 
non-commuting observables, increasing the uncertainty 
in one of the observables in order to reduce uncertainty 
in the other. In the case of quadrature squeezing, the two 
observables are amplitude and phase. Using typical no- 
tation, the squeezing in a state can be characterised by a 
parameter r, where the variances of the two quadratures 
are given by 



var(AT+) = 
var(A-) = e 2 



-2r 



(33) 
(34) 



For our purposes, the relevant property of this state is 
that it exhibits number squeezing. The expectation value 
and variance for the number operator N — a^a are given 



8 



by 

(N) = |a| 2 +sinh 2 r, (35) 
var(TV) = |q coshr — a*e~ t6 sinhr| 2 

+2 cosh 2 r sinh 2 r (36) 

where a is the coherent amplitude of the squeezed state. 

In order to stochastically simulate an atom laser 
sourced from a quadrature squeezed BEC, we require ini- 
tial noise that represents such a squeezed state, albeit the 
multimode rather than single mode version. Analogous 
to Eq. (|26|) . the Wigner distribution for a squeezed co- 
herent state is given by 

2 

W(xi,X2) = — exp \— 2 ({xi — a r ) 2 e~ 2r 

+ (x 2 - ai ) 2 e 2r )] . (37) 

This noise is applied to the initial fields at each of our 
grid points as before. For our simulations we chose a 
squeezing parameter r — In 2, which results in a num- 
ber variance of var(TV) = 7V/4. This reduction in the 
number variance by a factor of four over that of a coher- 
ent state corresponds to reducing the standard deviation, 
and hence the linewidth, by a factor of two. 

Figure [5] shows the linewidth of the atom laser in the 
long time limit over a range of BEC atom numbers. The 
results of stochastic simulations calculating the linewidth 
for a laser sourced from both squeezed and non-squeezed 
condensates are shown, along with the theoretical predic- 
tion for the linewidth based on Eqs. © and (f3"6")) . Figure 
|8] demonstrates that the linewidth limit scales as N 1 ^ as 
predicted, and that the reduction in number uncertainty 
due to number squeezing reduces the linewidth limit in 
agreement with Eq. (|36p . 

V. CONCLUSIONS 

We have quantified the lower bound for the linewidth 
of an atom laser due to the quantum diffusion of the BEC 
mode. We have simulated this process in one and two di- 
mensions, showing that the resulting linewidth limit is 
slightly lower than that estimated from a single mode 
BEC undergoing phase diffusion, but that the result 
scales in the same way as the zero-dimensional model. 

Just as the semiclassical model exhibits complicated 
spatial behavior when the physical timescales of the out- 
coupling and physical transport of the atoms from the 
outcoupling regime become comparable, so we might ex- 
pect significant spatial effects due to the quantum noise 
in the system as the coupling rate is increased. This rich 
behavior due to the interaction of the three timescales in 
the problem is essentially undesirable when considering 
an atom laser as an atom source, effectively resulting in 



an unpredictable atomic flux and energy spectrum. It 
will therefore be desirable to find ways to minimize this 
linewidth limit while remaining in the low outcoupling 
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Condensate number 



FIG. 8: Lower bound on linewidth as a function of condensate 
number for squeezed and non-squeezed condensates. Solid 
lines represent stochastic simulations; dashed lines the pre- 
diction according to Eqs. ([8]) and (|36[) . The upper pair of 
lines assume a non-squeezed condensate; the lower pair rep- 
resent a condensate with quadrature squeezing of magnitude 
r = 0.69. 



rate limit. This is a complicated trade-off: pumped sys- 
tems tend to be more stable and narrow when operating 
at high flux, and high flux implies either high atom den- 
sities (and therefore interactions), or high outcoupling 
rate, where the linewidth limit is higher. 

One obvious solution might be to attempt to reduce 
atomic interactions via some process such as a Feshbach 
resonance [32], but this introduces a further complica- 
tion, which is that earlier work has suggested that pump- 
ing schemes will tend to lead to modal instabilities ex- 
cept in the regimes of high interactions [U Gil ■ Further 
studies are required in order to determine the optimum 
way of finding a compromise between these competing 
processes. By the very nature of this compromise, these 
studies will require a quantum statistical model of the 
atom laser system with pumping that includes at least 
one spatial dimension. The model and results in this 
paper are an important step in that direction. 
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